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Teats  of  fit  for  the  Cauchy  distribution  based 
on  the  empirical  distribution  function 

by 
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Abstract 


Points  are  given  for  testing  'goodness-of-f it  to  the  Cauchy 
distribution,  with  unknown  location  and/or  scale  parameters.  The  tests  are 
based  on  the  empirical  distribution  function,  and  the  asymptotic 
points  round  off  work  begun  by  Darling  ( ]  95  5)  on  the  asymptotic 

theory  of  test  statistics.  Monte  Carlo  points  are  given  for 
finite  n  and  some  discussion  of  power  is  included. 


Key  words:  Goodness-of-f it  Tests 
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1.  INTRODUCTION. 

In  a  pioneering  paper.  Darling  (1955)  discussed  the  asymptotic 
theory  of  the  empirical  process  and  of  certain  goodness-of-f it  statistics 
based  on  this  process,  when  parameters  must  be  estimated  from  the  sample 
used  in  testing  fit.  The  estimated  parameters  were  location  and  scale 
parameters,  and  the  theory  was  illustrated  by  a  test  for  the  Cauchy 

2 

distribution.  The  statistics  discussed  were  the  Cramer-von  Mises  W  and  the 

2 

Anderson-Darling  A  ,  statistics  based  on  the  empirical  distribution 
function  (EDF)  of  the  given  sample. 

In  this  article  we  develop  the  tests  for  the  Cauchy  distribution, 

when  either  or  both  of  the  location  and  sclae  parameters  are  estimated  by 

efficient  estimators  given  below.  The  tests  are  set  out  in  Section  2. 

2  2 

Asymptotic  percentage  points  are  given  for  W  and  A  ,  and  also  for  the 
2 

EDF  statistic  U  introduced  by  Watson  (1961) ;  they  involve  calculating 

weights  in  sums  of  weighted  chi-square  variables.  This  is  done  by  techniques 

drawn  from  Darling  (1955)  and  the  details  are  given  in  Section  3.  For  finite 

samples,  points  for  the  three  statistics  have  been  found  from  Monte  Carlo  samples. 

Points  for  the  well-known  Kolmogorov  statistic  D,  and  for  the  related  V  were  found  at  tht 

same  time,  and  a  table  for  Case  3  is  given  for  reference;  the  asymptotic  theory 

used  for  the  Cramer-von  Mises  statistics  cannot  be  applied  to  these  statistics. 

2  2 

D  is  usually  not  as  powerful  as  W  or  A  ,  although  V  is  sometimes 

2 

competitive  with  U  .  A  brief  discussion  of  alternative  statistics,  and 
power,  is  given  in  Section  4. 
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2.  TESTS  FOR  THE  CAUCHY  DISTRIBUTION. 


Suppose  a  given  random  sample  is  X, ,X_,...,X  ,  with  order 

12  n 

statistics  X  <  X  <  ...  <  X.  .  .  The  test  discussed  is  a  test  of 
ll)  (2)  (n) 


Hq:  the  X-sample  comes  from  the  distribution 


F(x;o,B)  =  j  +  i  tan  1  ,  -  00  <  x  < 


with  density  function 


f(x;a,B)  = 


ir{l  +  (x-cO/B}" 


—  <*>  <  x  <  00 


We  can  distinguish  4  cases,  following  Stephens  (1974)1 


Case  0:  parameters  a  and  B  in  (1)  are  both  known. 


Case  1:  parameter  a  is  not  known,  B  is  known. 


Case  2:  parameter  a  is  known,  B  is  not  known. 


Case  3:  parameters  a,  B  are  both  unknown. 


In  Cases  1,  2  and  3  estimates  of  a,  8  are  obtained  from  the  formulas 

•  ~ 

a  =  E.  a.  X,..  and  8  =  I.  d,  X...  ,  where 
i  i  U)  ii  (i) 

sin ( 4tt {  i/ (n+1)  -  0.5}] 

Ci  n  tan { i/ (n  +  1)  -  0.5}] 


d^  =  8  sin[Ti{  i/ (n+1)  -  0.5}]cos  [ir{i/(n+l)  -  0.5}]/n  . 


Here,  and  in  later  formulas,  sums  run  for  i  from  1  to  n  . 
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These  are  not  maximum  likelihood  estimators  (MLEs)  but  are  asymptotically 
efficient,  an  important  requisite  for  the  asymptotic  theory  of  Section  3 
to  be  valid.  These  estimators  are  used  because  MLEs  are  known  to  be 
difficult  to  work  with;  the  likelihood  can  have  local  maxima  and  it  is 
sometimes  difficult  to  decide  on  the  global  maximum.  When  the  estimates 
are  obtained,  the  test  continues  with  the  following  steps: 

(1)  Calculate  z(.}  =  F(X(i);a,6),  replacing  a  and/or  0  by  estimates  where 
necessary ; 

(2)  Calculate  the  three  test  statistics  from 

W 2  =  ZjU  -  (2i-l)/(2n) }2  +  1/ (12n)  (5) 

U2  =  W2  -  n(z  -  0 . 5) 2  ,  where  z  =  ^z  j/n  (6) 

A2  =  -n  -n"1  Zi(2i-l){log(z(i))  +  log(l  -  z(n+1_i))}.  (7) 

Here  log  refers  to  natural  logarithm. 

(3)  Refer  the  value  of  the  statistic  used  to  Table  1,  for  the  appropriate 

Case:  is  rejected  at  significance  level  a  if  the  test  statistic 

exceeds  the  value  given  for  the  sample  size  n  and  the  desired  level  a  . 

The  present  Table  for  Case  0  is  a  more  accurate  update  of  a 

previously  published  table,  (Stephens,  1974,  1976)  although  the  changes 

2  2 

are  trivial  in  practice.  For  other  Cases,  the  distributions  of  W  ,  U  and 
2 

A  do  not  depend  on  the  true  a,  8  .  The  asymptotic  points  are  calculated 
from  the  theory  in  the  next  section,  and  points  for  finite  n  are  based  on 
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Monte  Carlo  studies  using  10,000  samples  for  each  n  .  It  can  be  seen  that 
for  this  very  heavy-tailed  distribution  ,  2nd  with  these  estimators  the 
points  vary  with  n  somewhat  surprisingly;  for  other.. distributions 
(see,  e.g.,  Stephens,  1974,  1977,  1979),  and  with  MLES,  they  converge 
much  more  rapidly  to  the  asymptotic  points. 


The  statistics  D  and  V  are  obtained  from  the  z..x 

(r) 

D+  =  max.  {  (i/n)  -  z  }  ;  D  =  max.  {z  -  ( i— 1 )  /n}  ; 
i  (r)  x  (r) 

D  =  max(D+,D  )  and  V  =  D+  +  D 


by 


Monte  Carlo  points  for  D/n  and  for  v/n  are  given  in 

2  2 

Table  2,  based  on  the  same  10,000  samples  as  for  the  statistics  W  ,  U 
2 

and  A 


J 
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3.  THEORY  OF  THE  TESTS. 

The  asymptotic  distribution  of  any  one  of  the  three  statistics 

is  that  of 


S  =  Ei  ui/Xi  ,  i  =  1,2,...  (8) 

2 

where  u^  are  independent  _x  variables  and  A^  are  weights.  The  weights 

are  found  from  the  now-classical  asymptotic  theory  of  the  empirical  process 

of  the  z- values .  Darling  (1955)  gave  this  theory  for  tests  for  absolutely 

2 

continuous  distributions  and  illustrated  it  for  W  with  the  Cauchy 
distribution,  although  details  of  how  the  A^  are  calculated  were  omitted 
except  for  Case  2.  We  now  complete  the  calculations,  following  the  steps 
and  notation  given  in  Stephens  (1976,  1977).  The  empirical  process,  for 
all  cases ,  becomes  asymptotically  a  Gaussian  process  Z(s),  with  E(Z(s))  =  0, 

Z (0 )  =  Z ( 1)  =  0,  and  with  the  covariance  p(s,t)  =  E(Z(s)Z(t))  varying  with 
the  Case.  In  Case  0  p(s,t)  =  pQ(s,t)  =  min  s,t  -  st  .  For  the  other  three  cases 
p(s,t)  takes  the  following  form: 

Case  1:  p(s,t)  =  pQ(s,t)  -  ^ (s) 4^ (t) 

Case  2:  p(s,t)  =  pQ(s,t)  -  <)>2(s)  <f>2(t) 

Case  3:  p(s,t)  =  pQ  ( s,t)  -  $  (s)  $  (t)  -  $2  (s )  $  (t)  with 

2 

^(s)  =  -  y/2  (sin  tts)/tt  and  <^2  (s )  =  (sin  2tiu)  / (  »£tt)  . 

These  results  for  Cases  1  and  2  were  given  by  Darling:  the  simple 
result  for  Case  3  follows  because  the  estimates  of  a  and  8  are  asymptotically 
independent  and  the  Fisher  information  matrix  is  diagonal  (see  Stephens,  1976, 
1977).  For  Cases  1  and  2  the  weights  A.  are  found  as  follows.  Firs*-  calculate 
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aj  =  j  ^1^  Si-n  ^ 


b_.  =  |  (}>2  (s)  sin  ffjs  ds  ,  j  =  1,2, _  . 


and  define 


S  (X)  =  1  +  X  Z 


l  - J-— -  ,  S  (A)  =  1  +  X  Z  - - 

>«1  1-X/Or  j  )  j=l  l-A/Orj  } 


Let  d^ (X)  be  the  Fredholm  determinant  associated  with 
2  2 

Pq  (s  ,t) :  dg(X)  =  H(X  -  ir  j  )  .  For  Case  1,  the  Fredholm  determinant 
3 

is  D  (X)  =  d  (X)  S  (X)  and  for  Case  2  it  is  D_  (X)  =  d  (X)  S,  (X)  .  The 
10a  2  0  b 

weights  for  these  Cases  axe  found  by  solving  D^CX)  =  0  for  Case  1  and 

D2 (X)  for  Case  2. 

Case  1.  It  is  easily  shown  that  a^  =  0  for  j  even,  and 
2  2 

=  8/{tt  j  ( j  -4)}  for  j  odd.  Setting  D^(X)  =0  gives  a  set  of 
♦  2  2 

solutions  X.  =  n  j  ,  j  =  2,4,6,...  ;  another  set  is  found  by  solving  S  (X)  =0 

1  cL 

2  .2 

(the  solutions  X^  =  tt  j  of  d^(X),  for  j  odd,  axe  not  solutions  of  D^(X)  =  0 

because  of  cancellation  with  the  denominators  in  S  (X)).  To  solve  S  (X)  =  0  it  is 

a  a 

00 

best  to  write  K  =  1/X  and  solve  S*(k)  =  1  +  Z  a2/{<  -  1/ (tt2  j 2 )  }  =  0; 

3  j=l  3 

a  solution  exists  in  each  interval  (1/(32tt2),  1/tt2)  ,  ( 1/ (5 2rr 2)  ,  1/(32tt2)), 

etc.  and  these  are  easily  found  numerically. 


Case  2 .  For  Case  2,  b,  =  0  except  b  =  l/2t  .  The  solutions  of  D  (X)  =  0 

D  *  * 
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2  2 

are  then  A^  =  IT  j  except  for  j 
remarked  on  by  Darling  as  losing  a 
of  $  . 


=  2.  This  rather  curious  result  is 
"degree  of  freedom"  by  the  estimation 


Case  3.  For  Case  3,  the  Fredholm  determinant  becomes 

D,  (X)  =  d„(X)S  (A)S,  (X),  (Stephens,  19  76)  and  setting  D  (A)  =0  gives  two  sets  of  A  ; 

3  0  a  b  3 

the  set  X*  of  S^tX)  =  0  already  found  as  part  of  the  solution  for 

Case  1,  and  the  set  A**=  l/n2 * * Sj2,  for  all  j  except  j  =  2  ,  found  for  Case  2. 

Cumulants  of  asymptotic  distributions.  The  cumulants  of  the  distributions 

can  be  found  by  direct  calculations  The  mean  for  Case  j  is 
1 

[  2  2 

It  =  1/6  -  I  <j>.(s)ds,  j  =  1,2  ;  the  values  are  y^  =  1/6  -  3/(4tt  )  =  0.0907 

0 

2 

(note  a  misprint  in  Darling,  1955,  Section  8a)  and  y2  =  1/6  -  1/(4tt  )  =  0.1413. 

2 

For  Case  3,  y^  =  1/6  -  1/IT  =  0.0653  .  Other  cumulants  may  be  calculated 

as  described  in  Stephens  (1976).  The  values  are  given  for  reference  in 
Table  3.  They  may  be  used  to  provide  checks  on  the  calculations  of  the  A.  , 
since  they  may  also  be  calculated  from  the  distributional  form 

S  =  E.u./X.  .  The  r-th  cumulant  is  =  2r  ^ (r-1) '.  £.  l/(A.)r  •  these 

ill  2  i  i 

converge  sufficiently  fast,  for  r  >  2  ,  to  give  excellent  checks  on  the  X^ 
values  by  matching  with  the  direct  calculations. 

When  the  X^  were  found,  for  the  different  cases,  Imhof's  (1961) 
method  was  used  to  give  the  percentage  points  for  S  .  The  points  were 
checked,  with  excellent  agreement,  by  fitting  Pearson  curves  to  the 
distribution,  using  the  first  four  cumulants.  The  slight  changes  from 
earlier  tables,  for  Case  0  points  in  Table  1  are  due  to  replacing 
Pearson  curve  fits  by  points  found  from  the  Imhof  method. 
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2  2 
Statistic  U  .  The  asymptotic  distribution  of  U  is  that  of 

1  1 

|  Z^(t)dt  where  Z^t)  =  Z(t)  -  j  Z(t)dt.  (Watson,  1961).  The  solutions 
0  0 

for  Ai  are  somewhat  more  complicated  in  principle  (see  Stephens,  1976) 

but  in  fact,  for  the  Cauchy  distribution,  they  work  out  easily;  details 

*  2  2 

will  be  omitted.  For  Case  1,  the  weights  are  the  set  A_.  =  4ir  j  ,  j  =  1,2,... 

*  *  *  it 

and  a  second  set  A.  which  are  identical  to  A.  except  that  *A  is  omitted. 

3  3  1 

For  Case  2,  the  weights  work  out  to  be  the  same  as  those  for  Case  1,  a 

2 

surprising  result  which  means  that  the  asymptotic  distribution  of  U  is 

the  same  in  both  Cases.  This  occurs  also  for  the  logistic  distribution; 

*  * 

see  Stephens  (1979).  For  Case  3  the  weights  are  two  sets  of  A.  .  The 
calculations  for  cumulants  give  the  values  in  Table  3. 

2  2 

Statistic  A  .  For  A  the  process  Q(t)  =  Z(t)/oj(t)  must  be  examined, 

where  w(t)  =  (t(l-t)}  .  The  details  parallel  those  given  in  Stephens  (1976) 
for  the  normal  distribution.  For  Case  0,  the  weights,  solutions  of  the 
corresponding  Fredholm  determinant  dQ (A) ,  are  A.  =  j (j  +  1) ,  j  =  1,2,...  . 

For  Case  1,  the  a_^  work  out  to  be  zero  for  j  even  and  must  be  found 

* 

numerically  for  j  odd  .  The  weights  A.  are  then  the  set  A  =  j  (i  +  1)  , 

3  3 

*  * 

j  =  2,4,6,...  and  a  second  set  A.  which  are  solutions  of  S  (A)  =  0  . 

3  a 

For  Case  2,  b  =0,  for  j  odd, and  the  weights  are  the  set  A*  =  j (i  +  1) , 

3  3 

*  * 

j  =  1,3,5, —  ,  and  th^  second  set  A.  ,  solutions  of  S^ (A)  =  0  .  For 

3  b 

Case  3  the  weights  are  the  two  sets  A**  for  Cases  1  and  2.  The 
calculations  for  cumulants  give  the  values  in  Table  3. 
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4.  FINAL  REMARKS. 

There  are  not  many  tests  available  for  testing  fit  to  the  Cauchy 
distribution.  In  this  article  we  have  given  points  for  EDF  tests,  tests 
which  are  consistent  and  unbiased  and  which  for  other  distributions,  are 
very  effective  in  terms  of  pcwer. 

Possible  alternative  tests  might  be  made  using  the  correlation 

coefficient  of  the  X^,  against  nt  ,  where  nr  is  the  expected  value 

of  the  ith  order  statistic  of  a  sample  of  size  n  from  (1),  with  Ct  =  0 

and  3=1.  The  values  of  nr  are  not  easily  obtained,  and  nr  might 

therefore  be  replaced  by  H  =F  ^(r;0,l)  with  r  =  i/(n+l)  .  H.  is  a 

i  1 

well  known  approximation  for  nr  for  most  distributions;  the  approximation 

is  less  good  in  the  tails,  and  of  course  for  the  Cauchy  distribution 

the  tails  will  be  important.  However,  is  easily  calculated  and 

tables  based  on  the  correlation  coefficient  between  X,.,  and  H.  ,  called 

U)  i 

R(X , H) ,  have  been  given  by  Stephens  (1986).  The  tables  are  for 
2 

Z(X,H)  =  n(l-R  (X,H) ) ,  and  are  given  for  complete  and  also  for  right- 
censored  samples.  other  possible  approaches  to  testing  fit  include 

tests  based  on  spacings  and  tests  based  on  the  empirical  characteristic 

function.  It  is  hoped  to  develop  such  tests  for  practical  use,  and  to 

include  them,  with  EDF  and  correlation  statistics,  in  an  extensive 

power  study.  Preliminary  work  suggests  that  EDF  statistics  are  much 

2 

better  than  correlation  statistics,  at  least,  with  U  and  V  best 
overall . 

The  author  expresses  thanks  to  the  Natural  Sciences  and  Engineering 
Research  Council  of  Canada,  and  to  the  U.S.  Office  of  Naval  Research,  for  support 


for  this  work. 
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Table  1 


Upper  tail  percentage  points  for 


U  and  A  . 


Significance  level  a  . 


n 

.25 

.15 

.10 

.05 

.025 

.01 

Case  1.  Statistic  Wz 

5 

.208 

.382 

.667 

1.26 

1.51 

1.61 

8 

.227 

.480 

.870 

1.68 

2.30 

2.55 

10 

.227 

.460 

.840 

1.80 

2.60 

3.10 

12 

.220 

.430 

.770 

1.76 

2.85 

3.65 

15 

.205 

.372 

.670 

1.59 

2-88 

4.23 

20 

.189 

.315 

.520 

1.25 

2.65 

4.80 

25 

.175 

.275 

.420 

.870 

2.10 

4.70 

30 

.166 

.250 

.360 

.710 

1.60 

4.10 

40 

.153 

.220 

.290 

.510 

1.50 

3.05 

50 

.145 

.200 

.260 

.400 

.70 

2.05 

100 

.130 

.170 

.210 

.270 

.35 

.60 

ao 

.115 

.146 

.173 

.216 

.260 

.319 

Case 

2.  Statistic  W* 

5 

.199 

.236 

.261 

.338 

.437 

.590 

8 

.211 

.273 

.321 

.389 

.463 

.564 

10 

.212 

.279 

.332 

.414 

.501 

.626 

12 

.212 

.281 

.337 

.433 

.525 

.661 

15 

.206 

.279 

.339 

.444 

.537 

.684 

20 

.199 

.273 

.333 

.442 

.547 

.698 

25 

.194 

.268 

.328 

.437 

.551 

.704 

30 

.189 

.265 

.326 

.435 

.553 

.708 

40 

.185 

.260 

.323 

.434 

.555 

.712 

50 

.183 

.258 

.321 

.433 

.557 

.714 

100 

.179 

.254 

.319 

.432 

.559 

.715 

OO 

.176 

.250 

.316 

.431 

.560 

.714 

Case 

3.  Statistic  Wl 

5 

.167 

.242 

.305 

.393 

.445 

.481 

8 

.192 

.315 

.441 

.703 

.940 

1.13 

10 

.197 

.331 

.481 

.833 

1.201 

1.571 

12 

.194 

.329 

.487 

.896 

1.391 

1.901 

15 

.185 

.317 

.472 

.904 

1.54 

2.33 

20 

.169 

.281 

.419 

.835 

1.63 

2.96 

25 

.154 

.253 

.366 

.726 

1.47 

3.08 

30 

.143 

.225 

.319 

.615 

1.25 

2.90 

40 

.126 

.195 

.263 

.460 

.850 

2.17 

50 

.117 

.175 

.235 

.381 

.642 

1.56 

60 

.1097 

.160 

.211 

.330 

.508 

1.07 

100 

.098 

.135 

.174 

.2378 

.331 

.544 

ao 

.080 

.108 

.130 

.170 

.212 

.270 

CASE  1.  STATISTIC  U2 


5 

8 

10 

12 

15 

20 

25 

30 

40 

50 

100 

oo 


.122  .173  .227 
.121  .185  .270 
.119  .180  .260 
.114  .172  .240 
.109  .158  .220 
.100  .141  .190 
.095  .128  .170 
.090  .121  .150 
.084  .110  .140 
.080  .104  .130 
.074  .095  .110 
.071  .088  .105 


.315  .387  .407 
.470  .600  .650 
.500  .720  .800 
.505  .780  .960 
.  480  .81  3  1.160 
.380  .780  1.340 
.280  .650  1.340 
.235  .480  1.230 
.195  .330  .970 
.170  .250  .600 
.145  .180  .250 
.133  .163  .204 


CASE  2.  STATISTIC  U2 


5 

8 

10 
1  2 
1  5 
20 
25 
30 
40 
50 
100 


.  1  20 

.  140 

.156 

.  1  22 

.  154 

.177 

.119 

.  149 

.  175 

.115 

.  144 

.  169 

.  109 

.137 

.  161 

.101 

.  126 

.  1  48 

.095 

.118 

.137 

.091 

.113 

.131 

.086 

.  105 

.  123 

.082 

.  102 

.117 

.076 

.096 

.111 

.011 

.088 

.105 

183 

.202 

.217 

221 

.280 

.358 

226 

.296 

.  400 

225 

.294 

.430 

210 

.276 

.403 

190 

.247 

.  355 

176 

.220 

.305 

166 

.202 

.270 

,154 

.  187 

.240 

,  148 

.180 

.230 

.  138 

.  169 

.210 

.  133 

.  163 

.204 

CASE  3.  STATISTIC  U2 


5 

.105 

.133 

8 

.107 

.  151 

10 

.  1  04 

.150 

12 

.  1  00 

.  1  44 

15 

.093 

.  1  32 

20 

.083 

.116 

25 

.075 

.101 

30 

.069 

.091 

40 

.062 

.079 

50 

.057 

.070 

60 

.054 

.066 

100 

.047 

.057 

CO 

.047 

.047 

.  160 

.202 

.226 

.198 

.293 

.386 

.  203 

.324 

.  461 

.200 

.339 

.504 

.  183 

.330 

.542 

.159 

.295 

.548 

.134 

.242 

.486 

.117 

.202 

.402 

.096 

.147 

.274 

.085 

.121 

.197 

.078 

.  104 

.  1  49 

.065 

.080 

.098 

.052 

.060 

.070 

.243 
.466 
.597 
.712 
.844 
.974 
.999 
.940 
.697 
.505 
.344 
.  154 
.081 


Case  1.  Statistic  A2 


s 

1.19 

2.22 

3.83 

8.00 

12.75 

17.980 

8 

1.33 

2.62 

4.7 

10.0 

17.4 

25.0 

10 

1.34 

2.52 

4.5 

10.6 

18.2 

29.0 

12 

1.31 

2.42 

4.1 

9.9 

18.8 

32.0 

15 

1.30 

2.15 

3.5 

8.2 

17.2 

31.2 

20 

1.17 

1.86 

2.8 

6.5 

14.4 

27.5 

25 

1.12 

1.68 

2.3 

4.7 

10.8 

23.0 

30 

1.08 

1.55 

2.1 

3.8 

8.2 

20.0 

40 

1.02 

1.38 

1.8 

2.9 

5.2 

15.5 

50 

.970 

1.29 

1.6 

2.4 

3.8 

10 

100 

.890 

1.16 

1.4 

1.8 

2.2 

3.5 

oo 

.834 

1.02 

1.219 

1.519 

1.812 

2.212 

Case  2.  Statistic  A2 

-■ 

5 

.974 

1.131 

1.239 

1.59 

2.08 

2.84 

8 

1.085 

1.360 

1.560 

1.88 

2.18 

2.55 

10 

1.110 

1.414 

1.653 

2.04 

2.38 

2.89 

12 

1.117 

1.443 

1.710 

2.14 

2.55 

3.15 

15 

1.117 

1.449 

1.728 

2.22 

2.65 

3.31 

20 

1.101 

1.444 

1.728 

2.24 

2.73 

3.44 

25 

1.083 

1.432 

1.727 

2.25 

2.77 

3.50 

30 

1.064 

1.422 

1.724 

2.25 

2.80 

3.53 

40 

1.051 

1.41 

1.723 

2.26 

2.82 

3.56 

50 

1.045 

1.405 

1.722 

2.27 

2.83 

3.59 

100 

1.038 

1.40 

1.718 

2.28 

2.86 

3.64 

oo 

1.034 

1.409 

1.716 

2.283 

2.872 

3.677 

Case  3.  Statistic  A2 

5 

.835 

1.14 

1.40 

1.77 

2.00 

2.16 

8 

.992 

1.52 

2.06 

3.20 

4.27 

5.24 

10 

1.04 

1.63 

2.27 

3.77 

5.58 

7.50 

12 

1.04 

1.65 

2.33 

4.14 

6.43 

9.51 

15 

1.02 

1.61 

2.28 

4.25 

7.20 

11.50 

20 

.975 

1.51 

2.13 

4.05 

7.58 

14.57 

25 

.914 

1.40 

1.94 

3.57 

6.91 

14.96 

30 

.875 

1.30 

1.76 

3.09 

5.86 

13.80 

40 

.812 

1.16 

1.53 

2.48 

4.23 

10.20 

50 

.774 

1.08 

1.41 

2.14 

3.37 

7.49 

60 

.743 

1.02 

1.30 

1.92 

2.76 

5.32 

100 

.689 

.927 

1.14 

1.52 

2.05 

3.30 

•O 

.615 

.780 

.949  1.225  1.52 

1.90 

Table  2 


Statistic  D 

n 

Upper 

tail  perc 

entage  points 

for  D 

and  V  , 

Case  3 . 

.25 

Significance  level  a  . 

.15  .10  .05 

.025 

.01 

10 

1.05 

1.22 

1.42 

1.75 

2.06 

2.37 

12 

1.00 

1.22 

1.42 

1.83 

2.22 

2.62 

20 

.946 

1.14 

1.32 

1.73 

2.25 

3.05 

30 

0.889 

1.05 

1.21 

1.54 

2.06 

2.98 

40 

0.850 

0.99  3 

1.12 

1.37 

1.77 

2.61 

50 

0.822 

0.949 

1.06 

1.28 

1.58 

2.29 

60 

0.802 

0.921 

1.02 

1.21 

1.42 

1.95 

100 

.755 

.755 

.925 

1.07 

1.23 

1.49 

Statistic ' 

n 

_V 

.25 

.15 

.10 

.05 

.025 

.01 

10 

1.30 

1.48 

1.65 

1.96 

2.27 

2.57 

12 

1.31 

1.48 

1.65 

2.01 

2.39 

2.79 

20 

1.24 

1.39 

1.53 

1.89 

2.36 

3.15 

30 

1.18 

1.30 

1.42 

1.69 

2.20 

3.09 

40 

1.15 

1.25 

1.34 

1.53 

1.91 

2.74 

50 

1.12 

1.21 

1.30 

1.46 

1.72 

2.40 

60 

1.10 

1.19 

1.26 

1.40 

1.47 

2.10 

100 

1.06 

1.14 

1.20 

1.30 

1.41 

1.64 

Table  3 


Cumulants 

of  asymptotic 

distributions 

10  y 

2  2 

0  X  10 

,3 
x  10 

4 

K„  x  ^0 

_ i__ - 

Case  0 : 

1.666 

2.222 

8.466 

50.79 

Case  1 : 

.9068 

.4052 

.5099 

.1015 

Case  2 : 

1.413 

2.094 

8.336 

.5060 

Case  3 : 

.6585 

.2769 

.  3799 

.8187 

Case  0 : 

.8333 

.2777 

.2645 

.  3968 

Cases  1,2: 

.5800 

.1495 

.1345 

.  1992 

Case  3: 

.  32  7 

.0211 

4.51  x  10"3 

1.58  x 

y 

o2 

*3 

k4 

Case  0 : 

1 

.5797 

1.043 

3.040 

Case  1 : 

.66  38 

.1872 

.1579 

.2137 

Case  2 : 

.8422 

.5249 

1.006 

3.003 

Case  3 : 

.5060 

.1324 

.  1211 

.  1768 

10 
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